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ABSTRACT 

A computational and experimental investigation has been 
undertaken to study the heat transfer characteristics of jet 
vanes used in thrust vector control systems. A computational 
model based on the PHOENICS code is being developed to predict 
the surface distribution of heat transfer coefficient over a jet 
vane immersed in a high temperature, high speed rocket exhaust. 
An experimental program is also being developed to verify the 
computational model. 
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JET VANE HEAT TOANSFER MODELLING 



INTRODUCTION 



jl£l_Yan£_lYC_Sy:5.i£ai5. 

Thrust vector control (TVC) systems offer means of tra- 
jectory control that are virtually independent of external 
forces. Such capability is frequently required for tactical 
missiles, as well as spacecraft launch vehicles, when the rel 
ative flow past external lifting surfaces is insufficient to 
generate the necessary control forces. This commonly occurs 
during low-speed operations, such as at launch or during ho- 
vering flight. High angle of attack flight may also lead to 
regimes where conventional lifting surfaces are inadequate. 

In addition, there are occasions when external steering de- 
vices are infeasible from a design point of view, such as for 
tube-launched devices. 

Several methods of TVC have been developed and applied to 
operational and experimental vehicles [1,2]. These include 
movable nozzles, internal fluid injection (secondary injec- 
tion), and mechanical jet deflection systems. Jet vane sys- 
tems fall in the latter category and they tend to be favored 
for volume-limited applications requiring relatively low actu 
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ation torques and rapid response. Among the various methods 
of TVC, jet vanes have been found to be the most effective in 
delivering large thrust deflection angles (up to about 30°), 
They may also be used with relative ease to generate roll 
torques . 

The application of jet vane TVC dates back to the rockets 
designed by Goddard, and has extended to the Redstone, Ser- 
geant, Talos, Pershing, and Algo II and III motors [1,2], as 
well as several intallations in smaller tactical rockets. 

Of course, there are disadvantages accompanying the 
selection of jet vanes for TVC purposes. These include thrust 
losses on the order of 3-5% with undeflected vanes [2]. In 
addition, the attainment of relatively high thrust deflection 
angles may lead to axial thrust losses of the same order of 
magnitude as the resulting side force. However, the chief 
problem associated with the use of jet vanes is the large 
thermal loading that they experience as they are required to 
operate in hot, high-speed, par t i c 1 e- 1 aden flows. This prob- 
lem leads to design limitations so that jet vanes are often 
restricted to shor t -dur at i on use in motors with low-tempera- 
ture non-metal i zed propellants. 

The aerodynamic (side-force producing) characteristics of 
jet vanes may be calculated with fair certainty on the basis 
of inviscid flow theory with suitable corrections for viscous 
effects [2,3]. On the other hand, difficulties that stem from 
the severity of the jet-vane thermal environment tend to lead 
to design practices that are based upon past experience and 
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cut-and-try methods. Over-design is therefore inevitable, 
with virtually no capacity for design optimization. In order 
to exploit the several advantages of jet vane TVC systems, 
therefore, it has become necessary to build reliable data 
bases and, to the extent possible, attain a fundamental under- 
standing of the heat transfer aspects of such systems. This 
has led to the work reported herein. 

In the design of a jet vane system, the integrity of the 
vanes themselves must be guaranteed over the specified work 
cycle. Although this is a serious challenge, designers must 
also also consider the behavior of the vanes and supporting 
structure during transient events. Upon motor ignition, 
localized jet vane temperatures may rise from their initial 
values to near-stagnation temperatures in fractions of a se- 
cond. Temperatures in the vane attachment device will rapidly 
follow this rise, and large thermal stresses may develop due 
to the proximity of a relatively cool supporting structure. 
These and other design aspects can only be addressed with pre- 
cision if there is a good understanding of the convective heat 
transfer process that gives rise to the energy transfer from 
the flowing gases to the the vane. Put another way, the 
application of complex computer codes for thermal conduction 
in the vane and supporting structure can only follow the 
specification of the convective boundary conditions. 

With respect to these conditions, the problem is even 
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more complex. The vane is immersed in a flow field that is, 
if generally described, compressible, turbulent, multi-com- 
ponent (and possibly mul t i -phase ) , three-dimensional, and 
unsteady, with variable properties and nonlinear and time- 
variant boundary conditions. These complexities, if taken more 
than one at a time, present a problem that is beyond the 
state-of-the-art for exact solution. In addition, the over- 
all flow field will contain intersecting and impinging shock 
waves that give rise to discontinuous events and further 
complication of the boundary conditions. The presence of 
various protuberances only serves to exacerbate these 
difficulties. 

To further define the problem, it may be of use to 
consider the levels of heat transfer that might be expected 
from the point of view of "simple" convection from a super- 
sonic flow to a cooled wall. The convective heat transfer 
coefficient (h) may be described in terms of the nond imens i on- 
al Stanton number (St) as follows; 

h = ( VCp) St 

where , V, and c^ are the density, velocity, and specific 
heat at constant pressure of the flowing gas. On the assump- 
tion that the gas behaves ideally, this expression may be 
written 

h = St (-^Y I g Igp) 

where P and T are local gas pressure and temperature and I 

s p 

is the specific impulse. Consider, for instance, the condi- 
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tions at the exit of a rocket nozzle flowing at Mach 3. With 
a ratio of specific heats assumed to be k = 1.2, we have, 
approximately: 

P 

h = St (0.25 g I ) 
o ^ 

where the subscript (^) refers to stagnation conditions. If a 
typical case is taken to be given by = 100 bar, = 2560K, 

and I = 250 s, then, in round numbers, 

sp 

h = 2000 St (kW/m^K) 

In rocket engine nozzles, a typical value for the Stanton 
number is about 0.002 [4]. Thus, with the result above, the 
expected value of the heat transfer coefficient is on the or- 
der of 4 kW/m K. Such a value should not be taken as conser- 
vative, since actual motor conditions may be more severe and 
local regions (e.g. stagnation points) may experience much 
higher values. By way of comparison, values of h such as 
these are well beyond those usually considered for convective 
heat transfer and, in fact, they are more typical of those 
realized in phase-change processes [5]. 

In short, there can be no argument that a complete and 
general analytical solution is too much to hope for. On the 
other hand, rapid advancements in the field of computational 
fluid mechanics give reason to expect that much progress can 
be made in predicting the major features of the flow field and 
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their effects on heat transfer to the vane. Even so, it will 
be necessary to provide considerable experimental backup in 
order to verify numerical predictions and to evaluate the 
effects of factors that are beyond present-day computational 
capabi 1 i t i es . 

This report describes the initial efforts in such a com- 
bined computat i onal -exper imental program. In the next section 
a brief review is given of some of the existing knowledge that 
may be brought to bear in the design of jet vane TVC systems. 
This is followed by a description of the code PHOENICS and 
discussion of some of the preliminary results that have been 
obtained through its use. The progress of the experimental 
program then will be summarized, again with some preliminary 
results. 
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BACKGROUND 



This section of the report has been prepared in order to 
provide an overview of some of the previous work that has been 
accomplished in the analysis of the heat transfer aspects of 
supersonic flows. The particular results that are discussed 
are meant to prove useful as "first-cut” design guidelines for 
the estimation of the thermal loads applicable to the jet vane 
configuration. With this goal in mind, an effort has been 
made to identify relationships that are backed by experimental 
results. Purely theoretical results will not receive atten- 
tion here. 

It is also intended that the reference list provided will 
serve as a starting point for those wishing to conduct a more- 
complete review of the available literature. Only the unclas- 
sified open literature is discussed here, with the result that 
some information of a practical nature may be omitted. On the 
other hand, it is felt that most proven results that are wide- 
ly applicable will have found their way into the public do- 
main. For those seeking information relative to particular 
applications, the Chemical Propulsion Information Agency has 
provided a useful bibliography [6]. 

In searching the literature (albeit in a rather cursory 
way), it has been observed that relevant articles seem to be 
grouped in time with large intervening gaps. A surge of in- 
formation appears in the mid-50s, in apparent response to 
stimuli provided by Sputnik and the emergence of supersonic 
flight as an accepted mode of transportation. The overall 
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impact of this work was to provide what was then seen as an 
adequate basis for the prediction of supersonic effects such 
as the aerodynamic heating of simple configurations. The com- 
plexity of the general thermal problem -- including, for in- 
stance, shock wave interactions, dissociation at high speeds, 
and complicated geometries resulting in three-dimensional 
flows -- was only accounted-f or by correction factors and 
overdesign. Since then, with a few noteworthy exceptions, the 
awesome nature of the complete problem has seemingly precluded 
major advances in the general understanding of heat transfer 
in supersonic flows. 

Since the mid-50s, much work has transpired in the under- 
standing of shock wave-boundary layer interactions. Secondly, 
the advent of space exploration has led to significant devel- 
opments in the state-of-the-art associated with hypersonic 
flight (with great emphasis on the space shuttle). With re- 
spect to experimental investigations, these advances have been 
greatly enhanced by the increasing use of facilities such as 
shock tunnels, gun tunnels, and Ludwieg tubes [7]. The third 
major factor associated with high-speed flows has been the 
explosive growth of the field of computational fluid mecha- 
nics. In this area, researchers are now able to predict flow 
effects that are most difficult, and sometimes impossible to 
measure with precision. From the point of view of this brief 
summary, neither of these latter two areas of advancement are 
particularly relevant. 

Insofar as thermal effects are concerned, there appears 
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to be a relative dearth of empirical design information that 
is directly and generally applicable to flows at supersonic 
Mach numbers below about six. For flows with intact boundary 
layers, the results of the mid-50s still seem to be appli- 
cable, with relatively few extensions since then. We hasten 
to repeat, however, that this conclusion is based on a view of 
the literature that is far from comprehensive. 

In the following discussion, initial emphasis is placed 
upon the thermal transport across an intact boundary layer 
separating a stationary flat plate from a supersonic free- 
stream flow. Following this, some departures from this more- 
or-less ideal situation are discussed, such as shock wave- 
boundary layer interactions and the effect of suface protu- 
berances . 

H£al_Iiiaii&l£ii-iQ-SucL£Ls.aiii$i-&ayLiida.Lit-Lait£iL-£I&ws.A 

The heat transfer to or from a bounding wall occurs on a 
microscopic scale at the wall. The rate of heat transfer is 
proportional to the temperature gradient at the wall, and in- 
creases or decreases according to variations in this gradient. 
For given temperatures of the wall and the external freestream 
(outside of the boundary layer), variations in the wall tem- 
perature gradient follow variations in the freestream velo- 
city. Increased speeds lead to thinning of the boundary 
layer, with the temperature change between freestream and wall 
occuring across the decreased thickness -- temperature gradi- 
ents and attendant heat transfer rates increase with free- 
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stream speed, and vice-versa. 

In high speed flows, these effects may be altered sig- 
nificantly by direct compression or viscous dissipation of 
energy by friction within the boundary layer. In such in- 
stances, fluid temperatures within the boundary layer may 
exceed freestream temperatures with attendant increases in 
heat transfer to a wall that is maintained below the free- 
stream temperature. In cases where the wall temperature is 
maintained above the freestream temperature, heat may be first 
transferred fr om the wall at low speeds and then la the wall 
as the flow speed is increased. 

In high speed flows, therefore, there is a certain wall 
temperature, above freestream temperature, at which there is 
no heat transfer to the wall. This temperature is called the 
"recovery" or "adiabatic wall" temperature. Whereas in low 
speed flows the freestream temperature is the reference tem- 
perature for heat transfer, the proper reference temperature 
in high speed flows is the recovery temperature. The recovery 
temperature, T^ , is defined in terms of the recover factor, r, 
as follows; 

T^ 

T = T, [1 + - 1) ] 

r 1 Tj 

where the subscript (j) is used to denote conditions at the 
boundary layer edge. Here, T^ is the maximum freestream tem- 
perature theoretically achievable (stagnation temperature): 
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so that 



Tj. = T^(l + r 

The calculation of heat transfer in high speed flows thus 
requires a determination of the recovery factor. The second 
important factor that di st inqui shes high speed flows from 
their incompressible counterparts is the fact that fluid pro- 
perties are, in general, dependent upon temperature. With the 
severity of temperature variations in such flows, and the pos- 
sibility of variations in density as well, the use of appro- 
priate fluid properties becomes a matter of great importance. 

As has been mentioned, a large amount of relevant infor- 
mation has been reported in the mid-50s timeframe. This is 
amply documented in a number of excellent review articles, 
among them Kaye [8], Eckert [9], and van Driest [10,11]. Ap- 
plications to rocket motors are discussed in Summerfield [12] 
and, more recently, Ziebland and Parkinson [4]. In the fol- 
lowing sections, we attempt to present the main features of 
these works. 

Lamina y Bounq^ry L&yqrSi For incompressible laminar 
boundary layer flow, application of the Reynolds analogy leads 
to a Stanton number given by: 

St = h/ CpV = 0.33/(Re^^^Pr^^^) 

It has been found that the form of this result can be carried- 
over into compressible flows E.iiay.id£Sl that the recovery tem- 
perature, as given above in terms of recovery factor and free- 
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stream Mach number and temperature, is used as the reference 

temperature for heat transfer. That is, the heat transfer 

rate from the wall, per -un i t -area , q’ ' , is given by 

q» « = h(T -T ) 

'1 ' w r 

In another form, 

q" = h(T„- Tj) tl - tj(r,Mj,T^/Tj)] 

where 






tzl 2 
r (=2--!-) 

- 1 
'Tl 



which illustrates that the modification for compressible flow 

depends upon both M, and T /T, . 

1 W 1 

For laminar flow, the recovery factor, and hence the 
recovery temperature, may be calculated with 1 % accuracy by 



for Prandtl numbers in the range 0.7 < Pr < 1.2 and Mach num- 
bers up to about 8 [8]. 

Use of these expressions further requires that fluid 
properties be evaluated at specified conditions. References 
[8] and [11] recommend the use of the freestream temperature 
for this evaluation, leading to the results shown in Fig. 1 
for air with a Prandtl number assumed constant at 0.75. Note 
that these results lead to a value of the "constant” coeffi- 
cient for the Stanton number that depends upon the Mach num- 
ber and the temperature ratio in compressible flow. The 
question of the correct reference temperature for the evalu- 
ation of fluid properties is further addressed by Eckert [9], 
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* 

in which this temperature is given the symbol T and is 
defined by: 




0.5 + 0.22 f^(r,Mj ^ 



where is defined above. 




Figure 1. Local Stanton number for a flat plate in air, 

Pr = 0.75 [8]. 

In summary, for high speed flows with laminar boundary 

layers the Stanton number may be estimated as follows: 

* *1/2 *2/3 

St = 0.33/(Re Pr ) 

where the star notation indicates that fluid properties are 

* 

to be calculated using T . From this, the film coefficient 
may be evaluated and hence the heat transfer rate per-unit- 
area, based upon the difference between freestream and 
recovery temperatures. 



IUCt?Ul£Dt. goundacy. Lay^CSA The details of turbulent 
boundary layer flows are poorly understood -- today, as they 
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were in the mid-50s. Turbulent flow heat transfer 



formulations, based upon the Reynolds analogy between skin 
friction and heat transfer are, of course encumbered by the 
same uncertainties. With the acceptance of their approximate 
nature, the adaptation of incompressible formulas to compres- 
sible flows follows the same basic philosophy as for laminar 
flows. That is, the reference temperature for heat transfer 

is the recovery temperature and the fluid properties must be 

* 

evaluated at T . 

For turbulent boundary layers, the recovery factor like- 
wise has been found to be mainly dependent upon the Prandtl 
number, as follows: 

D 1/3 
r = Pr 

The resulting recommended expressions for Stanton number are 

^ ^ 1/5 ^ 2/3 „ 

St = 0.0296/(Re Pr ), for 5x10 < Re < 10‘ 

and 



* *2/3 * 2.58 „ 

St = 0.185/[Pr (logj^gRe ) ], for Re > lO' 



These expressions assume a smooth plate surface. Modifica- 
tions for rough surfaces are necessary, but extensive system- 
atic information on rough surfaces exists only for the fric- 
tion coefficient in low speed flows. 

For both laminar and turbulent boundary layers, these 
expressions for high-speed flow are still recommended today 
for purposes of preliminary design [5]. 
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Qili£L_Eia£2.i££l_E£aiu.L£s._iii_sl£l_Yaa£_H£ai_liLaas.t£Ji 

The previous discussion has been related to the heat 
transfer associated with high speed flows with intact boundary 
layers in the absence of a pressure gradient. In many in- 
stances, at least locally, the boundary layer on a jet vane 
may be severely disrupted or even non-existent. A few refer- 
ences have been found that may be of use in predicting effects 
due to shock wave-boundary layer interactions, protuberances, 
and angle of attack. These are by no means the only complica- 
ting factors. Other details of the flow, which lie beyond the 
present coverage, include dissociation at high speeds and tem- 
peratures, surface ablation and other multi-phase multi-com- 
ponent effects, and boundary layer transition. 

Slifl2.!i_Way.£r.SQU.iidai.jL_LajL£]:_inl£r.ae.l.iflii5.A 

When the thermodynamic equilibrium of a substance is per- 
turbed, the effect of the disturbance is transmitted from 
molecule to molecule at the speed of sound. A new state of 
equilibrium is reached as the thermodynamic properties adjust 
to the requirements imposed by the disturbance. In an other- 
wise uniform flow, these requirements generally arise from 
changes in pressure that are typically caused by a deflection 
of the flow at a bounding surface or by the merging of two or 
more streams at different pressures. 

In subsonic flows, propogation of a disturbance at the 
speed of sound permits a gradual achievement of the new state 
of equilibrium in regions both upstream and downstream of the 
disturbance. In supersonic flows, on the other hand, fluid 
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elements cannot receive the disturbance "message" until they 
have reached or passed the source of the perturbation. While 
the necessary adjustments occur smoothly in a subsonic flow, 
they occur suddenly in supersonic flows -- instantaneously in 
the theoretical limit of an inviscid non-conducting gas. 

Disturbances demanding large pressure d,£C.£.£as.£S. in a su- 
personic flow are accomodated by sequences of small pressure 
decreases, and concomittant accelerations, that may be char- 
acterized as standing wave processes that are approximately 
isentropic. The region of adjustment is sometimes called an 
"expansion fan." On the other hand, when substantial pressure 
iaQ.L£as.£S, are required of a supersonic flow, partitioning of 
the necessary adjustment into a sequence of smaller pressure 
increases reveals that associated weak waves must coalesce 
into a single sudden event that is characterized by its name, 
the "shock wave". 

Substantial dissipation occurs within a shock wave and in 
the limit the compression and deceleration processes are math- 
ematically as well as physically discontinuous. For these 
reasons, the shock wave is not isentropic and its strength is 
directly related to the decrease in stagnation pressure of the 
gas passing through it. In a perfect gas, this decrease in 
stagnation pressure is easily related to the entropy increase 
across a shock. 

Flows that are at supersonic speeds relative to their 
surroundings cannot stay that way indefinitely. As such flows 
are brought to rest, either intentionally or naturally, the 
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attendant decelerations must be accompanied by shock wave pro- 
cesses. Thus shock waves are always "waiting in the wings" in 
supersonic flows, and the task of the designer of devices in- 
volving supersonic flows (such as rocket nozzles) is to ensure 
that shock processes will not be overly detrimental to perfor- 
mance. In supersonic diffusers, for instance, careful design 
can lead to gradual slowing of the flow through a series of 
relatively weak shocks. Thus the cumulative pressure increase 
is attained with less loss than that associated with a single 
shock process. 

The strength of shock waves (usually characterized by the 
concomittant pressure increase) depends upon the nature of 
their source. In the case of the deflection of an otherwise 
uniform supersonic flow by a rigid boundary, shock strengths 
increase with the deflection angle through the range of so- 
called "weak waves." These waves are oblique to the flow and, 
in inviscid flow, they emanate from the origin of relatively 
small deflections. Beyond a certain maximum deflection, de- 
pending on the upstream Mach number, shock waves detach from 
the deflection and become curved. In such "detached shock" 
cases, the flow downstream of the wave is extremely compli- 
cated. The curvature of the wave converts a one-dimensional 
upstream flow to a two-dimensional flow, with attendant non- 
uniformities in Mach number (including both subsonic and su- 
personic regions) and thermodynamic properties downstream of 
the shock. 

While only weak shock processes are observed when oblique 
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shocks are attached to a deflecting boundary, the detached 
shock exhibits the entire range of strengths from the normal 
shock, through strong and weak oblique shocks, to waves of 
vanishing strength far from the disturbance. In cases where 
boundary conditions require flow adjustments to surrounding 
pressures, the strengths of the shocks occuring in supersonic 
flows will depend upon the extent of the necessary pressure 
adjustment. In such free-jet flows, compressive shock pro- 
cesses are usually interspersed with alternating expansions, 
as in rocket exhaust plumes adjusting to changing atmospheric 
pressures . 

Finally, in this brief summary, it must be mentioned that 
shock and expansion waves seldom occur singly and without 
boundaries. When these waves arrive at rigid surfaces, they 
are changed in strength and, in most cases, reflected. Such 
occurences, involving shock intersections with boundaries, are 
the most common instances when inviscid approximations are 
likely to be unjustifiable. When waves intersect with one- 
another, there is a mutual modulation of their strengths. 

Flows downstream of such intersections are generally non- 
uniform and likewise subject to strong viscous effects. 

H£ai_lLaas.l£L_As.as.c.ls._fi£._Sliac.li_Wa5L£zEQU.Q!lar.y._LajL£Ji 
lal.£La£l.l£liaji. Boundary layers are happy in the presence of 
pressures that decrease in the direction of flow. When 
pressures increase, on the other hand, the pressure force 
supports the shear forces already active within the boundary 
layer. Such increasing pressures in the flow direction impose 
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what is called and "adverse pressure gradient," from the point 
of view of the boundary layer. Adverse pressure gradients 
promote transition from laminar to turbulent flow within the 
boundary layer, as well as its eventual demise -- separation 
from the boundary. Shock waves are pressure increase pro- 
cesses that are both sudden and extreme -- there are few 
events that are more "adverse" as far as the boundary layer is 
concerned . 

The boundary layer in supersonic flow is a transition 
region from external supersonic velocities to zero speed 
relative to the wall. Thus even in supersonic flow there are 
subsonic regions near the wall. Within these subsonic re- 
gions, it is possible to transmit pressure imbalances upstream 
so that the high pressures downstream of an impinging shock 
are felt in a region (rather than at a point) at the wall. 

The resulting influence on the boundary layer is more gradual 
than the shock process itself, and may extend for some dis- 
tance along a wall. Shock impingement on a flat plate may 
cause separation of a laminar boundary layer, followed by 
reattachemt with or without transition to turbulence. Tur- 
bulent boundary layers may or may not experience local 
separation, depending upon the strength of the impinging 
shock. Figure 2 illustrates these features. 

In the cases of more-complicated geometries, shock wave 
impingement may lead to permanent separation of the boundary 
layer. In any case, the effect is to greatly reduce the 
insulating properties of the boundary layer with local or 
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Figure 2. Examples of boundary layer behavior with an 
impinging oblique shock [13]. 



widely-spread decreases in resistance to heat transfer. Hung, 
et al . describe some shock wave-boundary layer effects 
observed in the space shuttle program [14]. Local heat 
transfer coefficents have been observed to be more than 100 
times those observed without shock impingement. In general, 
the boundary interference effects can be expected to lead to 
increased heating rates by factors of 15 to 20. 

In a jet vane system, shock wave patterns can be expected 
to be extremely complex. There are likely to be shock waves 
that emanate in the main nozzle and impinge on the vanes and 



20 



surrounding boundaries such as shrouds. These will, in turn, 
be reflected and secondary impingements may result. Even if 
the main nozzle and shroud do not generate significant 
shocks, the vanes will do so and these shocks may impinge on 
other vanes. Finally, the shock patterns will vary depending 
upon motor operating conditions and vane deflection angles. 

From a design point of view, however, it should be suf- 
ficient to identify the worst-case shock impingement problem. 
This may be quite difficult but, if it is correctly done, 
secondary and higher-order impingements should not lead to 
catastrophic events. For instance, a preliminary study might 
be conducted to identify the origin of the strongest possible 
shock during a jet vane excursion. Having done this, the 
point of nearest possible impingement could be predicted. 

Such a point would be a critical point for design purposes and 
the local shock strength there could be used to estimate the 
worst-case possibility with respect to convective heat trans- 
fer. 



Fortunately, there is an abundance of information in the 
literature that allows prediction of heat transfer amplifica- 
tion due to shock wave-boundary layer interaction [7, 14-16]. 
The most common method is to relate the ratio of heat trans- 
fer coefficients to the pressure ratio accompanying shock 
impingement : 



h ■ ^p '' 



In 



this expression, the subscript (. refers to the peak 
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values due to impingement while the subscript (^) refers to 
values calculated for the undisturbed boundary layer flow. 

The recommended value for the exponent n varies from 0.7 
to 0.85, depending upon the state of the boundary layer, with 
the higher value used for the turbulent case. The methods 
previously described are recommended for calculation of the 
undisturbed value of the heat transfer coefficient. It is 
important to note that the undisturbed pressure is not that 
behind the impinging shock, since the pressure rise at the 
surface is attenuated somewhat by the spread of the separated 
region. The references cited above will lead to appropriate 
values for the pressure ratio. 

The apparent validity of the correlation between heat 
transfer and pressure ratios in disturbed boundary layers is 
most fortunate since, in most cases, it is much easier to 
measure pressures than heat transfer coefficients. Further 
information relating to practical applications in this area, 
including heat transfer with ablation, may be found in [16] 
and [17]. 

F r PtUlpg rances . In many ways, the design configurations 
of jet vane systems give the impression that the vanes will 
appear as blunt bodies from a fluid mechanical point of view. 
In addition, practical designs often include structural mem- 
bers (fasteners, brackets, shrouds, etc.) that can appear as 
protuberances to oncoming supersonic flows. The nature of the 
interactions induced by such obstructions is basically one of 
a strong-shock upsetting a boundary layer. 
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Much work has been done in estimating the effect of pro- 
tuberances in boundary layers. Relevant information, as well 
as further references will be found in Whitaker and Ostowari 
[18] and Sedney and Kitchens [19], The latter work describes 
the geometry of the interaction regions for various obstacle 
shapes, and discusses the depen-dence of these regions on Mach 
number and Reynolds number. 

The heat transfer implications of protuberances in dis- 
cussed in Datis, et al. [20], and correlations are presented 
for estimating heat transfer amplifications. Highest ampli- 
fications occur usptream of protuberances, and can be as much 
as five-fold — on the order of stagnation point heating. 
Immediately downstream, regions of reduced heat transfer have 
been found with relatively mild increases observed far down- 
stream of the obstruction. 

Effect, of Angle Much of the preceding dis- 

cussion is only relevant to cases when the jet vanes are at 
relatively small angles of attack to the oncoming flow. In 
many instances, however, jet vane operation may occur at an- 
gles of attack up to about 30°. The effect of angle of attack 
has been considered by Motwani, et al. [21], in a series of 
tests of low aspect-ratio rectangular plates at various angles 
of yaw and pitch in subsonic flow under turbulent boundary 
layer conditions. 

The angle of yaw was found to have an insignificant 
effect upon heat transfer. The angle of attack effects were 
observed to be influenced by whether or not separation bubbles 
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were established at the leading edges of the plates. The p're- 
sence of leading edge separation, together with aspect-ratio 
(edge flow) effects, actually led to decreases in the average 
film coefficient below those expected for zero angle of at- 
tack. In addition, values decreased with angle of attack at 
the lower angles. At angles above 30°, dependency of heat 
transfer on angle of attack was no longer observed. It was 
noted, however, that the effects of increasing favorable pres- 
sure gradient with angle of attack may have been responsible 
for r el ami nar i zat i on of the flow with subsequent decreases in 
f i Im coef f icient . 

The tests described in [21] were complicated by a number 
of factors, and careful analysis would be necessary to deter- 
mine if these results are applicable to the jet vane problem. 
Expectations of decreased average heat transfer coefficient 
with angle of attack, based upon these results, are probably 
not warranted. In any case, high local heat transfer rates at 
leading edges are well-known occurences and, in the jet vane 
problem, leading-edge shocks are likely to have significant 
effects on the downstream boundary layer behavior. What can 
be said with certainty, is that further numerical and experi- 
mental work is necessary with respect to the effect of angle 
of attack on boundary layer flow and heat transfer. 
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COMPUTATIONAL MODELING 



A major objective of the work undertaken is to develop a 
computational tool which can be used to assess the heat transfer 
characteristics of jet vane thrust vector control systems. As 
described previously, the flow environment seen by the jet vane 
during operation is very severe; the impinging rocket exhaust 
flow is supersonic, very hot, probably with some solid particles, 
and possibly chemically reacting. Also, from a computational 
point of view, the geometry of the jet vane is highly complex. 
Finally, the flow is not planar but is more nearly conical. 

Only a few years ago the numerical analysis and calculation 
of such a flow field was completely impossible. In the past 
several years great strides have been made in the development of 
numerical techniques for computational fluid dynamics (CFD) . The 
impetus for this rapid advance in CFD has been the concurrent 
development of today's super computers. Indeed it is arguable 
that the needs of modern CFD codes has been a major cause for the 
development of super computers. Whatever the underlying causes, 
we have reached the point in the development and implementation 
of CFD that an engineering analysis of the jet vane heat transfer 
problem is a real possibility. 

These rapid developments in CFD are well documented in the 
recent books by Peyret and Taylor [22], Patankar [23], and 
Jaluria and Torrence [24]. The problems associated with the 
numerical analysis of compressible flows are covered in the 
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survey article by MacCormack and Lomax [25] while the compilation 
of papers edited by Book [26] provides a great deal of insight 
into the very specialized vectorizing techniques necessary to 
efficiently implement CFD on modern super computers. 

The state of maturity of CFD is evidenced by the fact that 
there are now commercially available large scale general purpose 
computer codes for analyzing complex flow situations. Examples 
of these include: FLUENT and its variations from Creare Inc., in 
Hanover NH; FIDAP from Boeing Computer Services in Seattle, WA; 
and PHOENICS from CHAM of North America in Huntsville, AL. Since 
the task at hand was to develop an engineering analysis and 
design tool for assessing the heat transfer characteristics of 
jet vane thrust vector control systems and not to undertake an 
extended program of CFD code development, it was decided to use 
one of the commercially available codes and adapt it to the jet 
vane geometry and flow conditions. 

Because of prior favorable experience at the Postgraduate 
School with earlier versions of CHAM codes it was decided to use 
PHOENICS to implement the computational modeling. The 
availability of Mr. Amiram Leitner, of the Isreali Ministry of 
Defense, who was spending a sabbatical year at the Postgraduate 
School added a note of serendipity to the choice of PHOENICS. 

Mr. Leitner is the leader of a software development group for the 
Isreali Defense Research Lab, and is an expert in the use of 
PHOENICS for the analysis of missile related flow problems. Mr. 
Leitner's expertise was invaluable in the computational modeling 
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effort of the past year. 
Description of the PHOENICS Code 



As described in the instruction manual: "PHOENICS is a 
general-purpose computer code for the simulation of fluid-flow, 
heat -t r ans f er , mass-transfer , and combustion processes. It 
handles one- or two-phase flows, in one, two or three dimensions; 
and it can be applied to numerous problems arising in 
engineering, scientific research and the environment." The word 
PHOENICS is an acronym for E.arabolic, Hyperbolic, Qr E.lliptic 
Numerical Lntergration Qode S.eries. (The reference to parabolic, 
hyperbolic or elliptic is to the partial differential equation 
class to which the Navier Stokes equations may belong depending 
on the magnitude of the Reynolds number and Mach number.) 

Built into PHOENICS are the equations , and the 
corresponding computer coding sequences, which express the 
conservation laws of; mass of each phase present; the three 
components of momentum for each phase; the thermal energy of each 
phase; the mass concentration of each chemical species; and the 
fluxes of thermal radiation. Also included are turbulence models 
which account for the conservation of the energy and length scale 
of the turbulence. The equations contain individual terms which 
express the processes of convection, diffusion, generation and 
dissipation, with appropriate account taken of the interaction of 
turbulence with laminar flow phenomena. The conservation 
equations are applied to a large number of contiguous sub-domains 
or cells, into which the region of study is artificially devided. 



27 



The PHOENICS computer code is modular, with separate modules for 
input, output, computation, restart, and other computer "house 
keeping” chores. A more complete description of the structure of 
PHOENICS and its use can be found in the documentation for 
PHOENICS [27]. 

Use of PHOENICS for the NWC Jet Vane C onfiguration 

The geometry of the NWC jet vane configuration is shown in 
figure 3. As is obvious the geometry of the vane is three 
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Figure 3. NWC Jet Vane Configuration 



dimensional and quite complex. The modeling of such a problem is 
even more difficult when it is recalled that the impinging flow 
is conical and the vane operates at angle of attack. To gain 
initial familiarity with the use of the PHOENICS code and at the 
same time gain useful engineering information about the surface 
heat transfer characteristics of the vane initial two- 
dimensional approximations of the vane geometry were constructed. 
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Two approximations were used, one with a sharp leading edge and 
one which attempted to approximate the actual rounded leading 
edge of the vane. The dimensions of these two approximations are 
illustrated in figure 4, 

All Dimensions in Millimeters 




Figure 4. NWC Jet Vane Approximations 

The major difficulty in numerically modeling such a complex 
flow situation is in constructing an appropriate finite 
difference grid to capture and resolve some of the complicated 
physical phenomena such as the viscous boundary layer over the 
surface, the leading edge shock wave, or the expansion waves from 
the shoulder and the base of the vane. The grid construction and 
the other details of the implementation of the PHOENICS code for 
the NWC jet vane approximations are completely covered in the 
report by Leitner [28] and the Master of Science thesis by 
Yukselen [29]. What will be presented here are some of the 
results of the PHOENICS calculations taken from these two 
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documents . 



Results of Calculations 

Four cases were calculated: sharp edge vane and blunt edge 
vane; laminar boundary layer and turbulent boundary layer on each 
vane. For all cases the free stream Mach number was 3.2. The 
results presented in figures 5 through 8, are given as streeunwise 
variations of local surface shear stress coefficient (C^), or 
local Stanton number (S^). Figure 5 gives the distribution of 
the shear stress coefficient over both the sharp edge vane and 
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Figure 5. for a Laminar Boundary Layer 
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the blunt vane with a laminar boundary layer. As would be 
expected the shear stress shows a maximum at the start of the 
vane and falls rapidly downstream with another jump at the 
shoulder of the vane. Over most of the wedge portion of both 
vanes the value of is independent of the geometry of the 
leading edge. The blunt vane shows a higher initial and a 
higher at the shoulder than for the sharp edge vane. Figure 
6, shows the distribution of Cf for the same two vanes with a 
turbulent boundary layer. The shear stress coefficient for the 
blunt vane starts much higher than the sharp edge vane but drops 
rapidly to a value significantly lower than the sharp edge vane 
over most of the vane surface. Figures 7 and 8 give the 
distribution of the local Stanton number for both vanes for both 
laminar and turbulent boundary layers. In both figures the 
trends are the same as for the distributions. From these 
results it appears that the influence of the leading edge 
geometry is much more pronounced for the turbulent boundary layer 
than for the laminar boundary. At the high speeds and the length 
scales involved in these cases it is expected that the boundary 
layer would be turbulent. 

The cross over that is observed to result due to the leading 
edge geometry with the turbulent boundary layer, is not fully 
understood at this time but numerical anomolies cannot be ruled 
out. The shape of the leading edge, of course, has a great 
effect on the rate of growth of the boundary layer and on the 
shape and position of the initial shock wave. It is felt that 
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the influence of the leading edge geometry on the shock 
probably the more significant phenomenon in causing the 
cross over. 



wave is 
observed 
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Figure 6. for a Turbulent Boundary Layer 
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Figure 7. for a Laminar Boundary Layer 
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Figure 8. for a Turbulent Boundary Layer 
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EXPERIMENTAL PROGRAM 



A program of experimental investigation has been begun 
with two main goals in mind. These are (1) to develop a capa- 
bility to experimentally evaluate the accuracy of the heat 
transfer predictions evolving from the CFD analyses and (2) to 
provide a means for estimating the effects of complicating 
factors that are beyond the state-of-the-art in numerical mo- 
deling of jet vane flow fields. 

Two flow facilities are available at the Naval Postgradu- 
ate School to support these tests. The low-speed low-turbu- 
lence wind tunnel provides flow at speeds up to 100 m/s in a 
test section 1.8-m long, 0.7-m high, and 0.5-m wide. A small 
blow-down supersonic wind tunnel is also in place, driven by a 
2.4-MPa (350 psia) supply system. The test section of this 
facility is approximately 10-cm square with nominal test-sec- 
tion Mach numbers of 1.5, 2.0, 2.8, and 4.0 attainable using 
interchangeable nozzle-blocks. 

H£ai_lLans.l£L_M£as.u.ii£ai£nis._in_Hig.lizS[i££d[_EIflWS. 

The experimental deduction of convective heat transfer 
characteristics necessarily involves the measurement of two 
basic quantities: the driving temperature difference and the 

resulting rate of heat flow to or from the wall. Under steady 
state conditions, these may be measured by noting the heat 
flux needed to maintain the wall at a specified constant tem- 
perature in the presence of the external flow. This classical 
method has been extensively used for low speed flows where 
large test surfaces can be evaluated with relatively easy 
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access for the necessary instrumentation and control [5,30], 

In high-speed flows, where models tend to be relatively 
small and inaccessible, a common method for heat transfer 
experimentation is the so-called transient method. In this 
method, models are designed with thermal characteristics such 
that internal resistance to heat conduction is negligible with 
respect to that of the convective film. In such a situation 
the model temperature distribution may be assumed to be rela- 
tively uniform and the major temperature gradient would be 
that obtained across the convective boundary layer. The re- 
sulting analytical method is one of ’’lumped heat capacitance” 
with the model taken to be a uniform temperature sink (or 
source) with respect to the convective boundary condition. 

From an energy balance for such a case, the rate of heat flow 
across the boundary is 



a’ = hA(T -T ) = - me 

w r dt 

where A is the convective surface area and m and c are the 
mass and heat capacity of the model. 

Solution of this simple differential equation leads to 
the prediction of a wall temperature response given by 
T - T 
i r 

where T. is the initial wall temperature. The time constant 
for the event is given by 



1/t 



c 



= 

me 
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This time constant may be further characterized through the 
use of the non-dimensional Biot and Fourier numbers: 

Bi = hv/kA, Fo = a (A/v)^ t 

where v is the model volume and a is its thermal diffusivity. 
The elapsed time may now be written in terms of time constants 
as : 

t/t^ = Bi X Fo 

The appearance of the Biot number in these expressions is par- 
ticularly informative since the value of this number expresses 
the resistance to conduction in the model relative to the 
convective resistance of the external flow. According to the 
premise of lumped thermal capacity, the value of the Biot 
number must be small -- suggested limiting values are on the 
order of Bi 1 0.1 [5, p. 101]. 

With experimental procedures suitably designed, the 
measurement of the model temperature variation with time leads 
to the determination of the time constant for the response and 
thence to the film coefficent, h. A number of practical dif- 
ficulties arise, including attainment of a desired initial 
model temperature and the maintenance of this temperature dur- 
ing wind-tunnel start-up transients. If the rate of heat flow 
along the model is not negligible, the model temperature may 
become nonuniform. Errors due to this source may be minimized 
by short measurement periods (in which the change in model 
temperature does not give rise to significant streamwise 
temperature gradients) as well as the maintenance of small 
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Biot numbers . 



In compressible flow, the transient method requires the 
calculation or measurement of the recovery temperature. If 
the run duration is sufficiently long, this equilibrium tem- 
perature may be extrapolated in time (providing that stream- 
wise conduction in the model is not a factor). Measurement of 
local model surface temperatures may also be a problem and, in 
such cases, surface temperatures are often deduced by means of 
thermocouples imbedded within the model -- the so-called 
inverse conduction problem. These and other factors are dis- 
cussed in Bleakney [31] and Burbank and Hodge [32]. Appli- 
cations using imbedded thermocouples are developed and refined 
in Williams and Curry [33], Mehta [34], and Raynaud [35]. 
Et£LimiiiaLji_Exii£Liai£nlaL_WaLli^ 

As noted above, the experimental program has two main 
goals. To verify (or refute) the results of the numerical 
calculations, tests were designed to measure heat transfer 
coefficients in supersonic flows past standard shapes. For 
the first tests the nominal test-section Mach number was 
selected to be 2.8 and the model configuration was a double- 
wedge (diamond) . 

MQdg,l - Design 1 in order to use the transient method, 
model design necessitates a configuration leading to the 
smallest possible Biot number. To do this, it becomes neces- 
sary to estimate the film coefficient for the nominal test 
cond i t i on . 

For a double-wedge with 6-degree half-angle in Mach 2.8 
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flow, the freestreem Mach number downstream of the attached 
oblique shock is about 2.52. With nominal stagnation condi- 
tions of 0.7 MPa (100 psia) and 290 K (520 R), corresponding 
freestream static conditions are = 38.7 kPa (5.6 psia) and 
= 127 K ( 229 R). The freestream velocity may be calculated 
to be 570 m/s (1870 ft/s). For standard air properties [36], 
the kinematic viscosity at these freestream conditions is 

found to be 8.3x10 ^ m^/s (8.9x10 ^ ft^/s), with a correspon- 

7 -1 

ding Reynolds number per-uni t-length of 2.1x10 ft . For 

models of the size contemplated for the small supersonic wind 

0 

tunnel, this leads to Reynolds numbers on the order of 4x10 
-- marginal with respect to the possibility of laminar flow on 
the model surface. 

Following the methods described earlier for compressible 
flow, the following values may be calculated for laminar flow 
at a point 5 cm along the wedge with a wall temperature of 
290 K: 

r = 0.85 

T = 211 K 
r 

T* = 240 K 

Pr* = 0.72 

♦ 6 
Re = 3x10 

St* = 2.4x10"^ 
h = 227 W/m'^K 

The corresponding value for film coefficient if turbulent flow 
is assumed is 1780 W/m K. This latter value is assumed to be 
the most likely in actual testing and, in any event, it fixes 
the worst case in terms of minimum Biot number. 
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Copper [k = 400 W/mK (230 Btu/hr-f t-R) ] was chosen as the 
model material. For this configuration, shown in Fig. 9, the 





_ 3 

volume-to-ar ea ratio is v/A = 4x10 m and, using the turbu- 
lent value estimated above for the film coefficient, the ex- 
pected nominal Biot number is about 0.018 -- well below the 
desired maximum. 
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Me.a5.iiLeai£uls._U5.ia&_lalLazL£d_Ili£i:j[iia£Laiiliii^ 

Application of the transient method for heat transfer 
experimentation requires the measurement of surface tempera- 
tures, preferably at every point of interest on the model. 

This presents several significant problems; (1) location of 
thermocouples on a surface is often difficult, and the wiring 
necessary may lead to interference with the thermal and flow 
fields, (2) response times of thermocouple temperature sensing 
systems may be so slow as to distort the recorded transient 
response, (3) prior to testing, investigators may not know 
what points are liable to be of particular interest -- good 
experimental insight and probably luck are necessary to pick 
up anomolous behaviors, and (4) pointwise measurements do not 
indicate special trends on the model surface -- observed 
effects may be due to extreme gradients not detected by a 
given thermocouple array. 

With recent IR sensing and recording technology, these 
difficulties are virtually eliminated. In addition, the 
second goal of the experimental investigation -- to test com- 
plex jet vane geometries in various orientations -- is especi- 
ally enhanced by IR thermography. 

The thermal imaging system that has been employed in the 
experiments described below is the Probeye Thermal Video Sys- 
tem, Series 4300 [37], This system provides real-time temper- 
ature maps of objects through 16 color' bands displayed on a 
color monitor and recorded on video tape. The temperature 
range of the system is from 233K to 553K (-40xF to 536xF) and 
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the 16 color bands may be selected to span all or portions of 
this range. Emissivity adjustments can be made from 0.0 to 

1.0 in increments of 0.01, in the spectral wavelength band of 

2.0 to 5.6 fm. The unit includes a movable cursor which can 
be placed at any location in the field of view, giving digital 
readout of the temperature at that location. A real-time 
clock gives readouts precise to 0.1 seconds. 

The use of this imaging system requires employment of an 
optical window transparent to IR radiation in the wavelengths 
to be sensed. For this purpose, a six-inch diameter window, 2 
mm in thickness and fabricated of magnesium fluoride has been 
used in experiments in the subsonic wind tunnel. This materi- 
al, tradename IRTRAN 1 (Eastman Kodak Company), allows over 
959$ transmittance in the 2.0 to 5.8 fm range [38]. The window 
is capable of withstanding differential pressures of up to one 
atmosphere and ambient temperatures, with a safety factor of 
four . 

Before the system for use in the supersonic wind tunnel 
could be designed, fabricated, and installed, this facility 
was deactivated pending safety tests. At the time of this 
writing, these tests have not yet been completed. According- 
ly, the subsonic tunnel was selected in order to continue the 
work and to gain familiarity with the test procedures and 
equipment, and to develop validated methods for data acquisi- 
and reduction. The details of these tests are recorded in 
Spence [39 ] . 
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The model designed for the supersonic tests was fabri- 
cated without heating provisions. Because of the accessibil- 

I , 

ity of the model when mounted in the subsonic tunnel, it was || 

I 

possible to heat the model using a thermal blanket device that 
was manually removed at the instant prior to starting the 
tunnel . 

a£S.uIis._lat_ili£_Dau.tl£zW£d£e.j. 

Following extensive testing of the details of the experi- | 
mental setup and calibration of the sensing instrumentation, 

I 

production runs were conducted over the range of available 
wind tunnel speeds. The transient method was used to reduce 
the data from the surface thermocouple (Fig. 9) as well as | 

from the IR indications given with the cursor overlayed on the ^ 
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Figure 10. Temperature-t ime responses meaured by thermocouple 
and IR sensor; 306 ft/s freestream speed. 
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thermocouple location. Figure 10 shows typical time histories 
of these temperatures, indicating the agreement of the tem- 
perature decay rates with expected exponential form. In the 
table below, results are shown for both sets of measurements. 
The experimental uncertainty for the two representative speeds 
shown was on the order of 2 Btu/hr-ft -R [39], indicating re- 
markably good agreement between the two methods -- virtually 



Tunnel Film Coefficient, h 

Speed 

ft/sec Thermocouple IR System 

200 27.5 - 33.2 27.7 - 33.5 

306 41.0 - 42.7 39.9 - 41.8 



identical within experimental uncertainty. The theoretical 
flat plate estimates for the film coefficient are 22 and 28 
Btu/hr-ft^-R for 200 and 306 ft/s, respectively, on the as- 
sumption of laminar flow. In both cases, this assumption is 
marginal, and the results indicate that turbulent or transi- 
tional flow probably existed at the location of the measure- 
ment. (This conclusion is further supported by the increased 
deviation at the higher speed.) 

These tests served as a valuable learning experience. In 
hindsight, after encountering and working-out many experimen- 
tal bugs, it seems probable that the planned initial testing 
in the supersonic wind tunnel was i 1 1 -conce i ved . In addition, 
these results validate the application of the transient method 
using IR thermography. 
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It is important to note that the measurements given above 
correspond to conditions at a particular point on the wedge -- 
the location of the thermocouple in this case. Further test- 
ing, with the cursor positioned at other points, would allow 
the development of a complete heat transfer map of the wedge 
surface. However, viewing of the thermal behavior of the 
entire wedge surface showed that the transient isotherm vari- 
ation was essentially uniform over a large portion of the 
surface (with the exception of narrow bands in the three-di- 
mensional boundary layer regions near the sidewalls). Thus it 
may be deduced that the values determined above, at the ther- 
mocouple location, are representative of conditions over a 
major portion of the wedge surface. The ability to make such 
observations, which might otherwise require the use of hun- 
dreds of simultaneous point measurements, is a significant 
virtue of the IR thermography method. 

I R Thermography of a Jet Vane. The advantages of the IR 
thermography method are especially dramatic when it is used to 
view the thermal behavior of complicated configurations lead- 
ing to flow effects that are well beyond exact analytical 
capabilities. Such is the case for the jet vane, and some 
preliminary testing was done in order to demonstrate the 
method . 

The vane tested, shown mounted in the tunnel in Fig. 11, 
was constructed of copper-impregnated tungsten with an 

-4 

estimated vo 1 ume- to-sur f ace area ratio of about 1.5x10 ft. 
Thus, even though an exact value could not be calculated, the 
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Figure 11. Photograph of jet vane model mounted in subsonic 
wi nd tunnel . 

Biot number was well within the acceptable range. This was 

verified by the exponential temperature responses obtained in 

the tests. At 200 ft/s, the midspan heat transfer coefficient 

2 

was measured to be 41.5 Btu/hr-ft -R, while a measurement at 
the vane/baseplate intersection indicated a 32 % increase over 
this value. 

By simply viewing the time history of the thermal re- 
sponse of the vane, it was possible to note the expected high 
heat transfer rates at the vane leading and trailing edges. 
Other areas of high heat transfer were noted just downstream 
of the ridge of the vane, particularly near the tip. In other 
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areas, the heat transfer appeared to be more-or-less uniform. 
Again, these tests proved the effectiveness of the method and 
indicated many possibilities for future research using this 
application of IR thermography. Further details regarding 
these tests, including color prints of the IR images, may be 
found in Spence [39]. 
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CONCLUSION 



Preliminary work has been carried out to develop the 
capability to model the heat transfer processes encountered by 
jet vane thrust vector control systems. Computational modeling 
has begun using the PHOENICS code. The initial configurations 
investigated were two-dimensional approximations to the NWC jet 
vane geometry. Two separate leading edge geometries have been 
investigated: a sharp wedge leading edge; and a rounded nose 
wedge. The results from these calculations indicate that the 
leading edge geometry has a very strong effect on the skin 
friction and surface heat transfer over the whole body. Work is 
continuing to investigate the effects of the leading edge 
geometry on the surface skin friction and heat transfer. It is 
believed that these effects stem from the complex shock wave 
boundary layer interaction near the leading edge. To carry out 
this investigation the latest version of the PHOENICS code 
(PHOENICS 84) is being procured. This version of PHOENICS will 
allow more precise modeling of various leading edge geometries as 
well as the modeling of three dimensional flow configurations. 

In support of the computational modeling effort an 
experimental program has also been initiated. The purpose of 
this effort is to obtain confirmation of the computational 
results. The initial experimental work has been carried out in 
the subsonic wind tunnel and it is intended to continue the 
experimental work in the NPS supersonic wind tunnel. 
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